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ABSTRACT 



The objective of this study is to model and simulate using an Analog 
computer the performance of a relatively complex engineering system. The 
results of the simulation are validated by a comparison of the simulation 
with the measured performance of the prototype system. 

The prototype system consists of a single shaft gas turbine engine 
driving a three-phase alternator. The alternator load is a variable speed 
three phase wound rotor induction motor which drives a ship bow thruster 
propeller. In order to provide sufficient data to validate the simulation, 
the alternator speed, current and terminal voltage, KW load, and induction 
motor speed are recorded continuously as the motor is switched through five 
speeds. The recorded transients in these variables are then used to validate 
the model and simulation. 

The system is segregated into five sections and analytical models are 
determined for each component. The gas turbine is modeled using a transfer 
function derived from the engine performance curves. The speed governing 
system is modeled using transfer functions for the governor components as 
supplied by Woodward Governor Company. The alternator model is a result of 
the work of Professor Paynter and is taken from notes from his subject 
"Dynamics and Control of Rotating Machine Systems." The induction motor 
model is derived from a simplified equivalent circuit. The propeller is 
modeled as a torque proportional to the square of the speed. The simula- 
tion is simplified considerably by using a single-phase model to represent 
the electrical machinery. 

The correlation between the simulation and the prototype performance 
is excellent. The models used are relatively simple and yet can provide 
close estimates of the system performance and can demonstrate how perform- 
ance varies as the system parameters are changed. 
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1. INTRODUCTION 



Simulation ha.s become one of the most important tools available 
to engineers in the study and design of complex physical systems. The 
ability to predict a system 1 s response to various inputs is useful in 
many ways. Actions which can have disastrous effects on a system can 
be avoided if accurate predictions of the system response are avail- 
able. Simulation can also aid in the realization of high levels of 
system performance by aiding in control system design. Studies of the 
effects of varying system parameters can be performed with a simulation 
so that an optimum or near optimum prototype can be built without a 
costly trial and error process. 

The first stage in simulating a complex physical system is to i- 
dentify the individual components. An attempt is then made to formulate 
a model for each component of the system so that the interconnected mod- 
els will predict the performance of the system within a specified level 
of accuracy. Simulation is in effect a four-stage process. In the 
first stage the system is reduced to individual components and "black 
boxes’ 1 are used to represent these components. In the second stage, 
the components are interconnected so as to reflect the structure of the 
prototype system. The third stage usually requires the use of a com- 
puting machine to determine the time history of the variables as the 
components interact with each other. 

This simulation procedure can be carried out almost completely us- 
ing the concept of Bond Graphs (1), (2). Although bond graph techniques 
can be used to develop a set of state equations representing the entire 
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system, for this simulation bond graphs will be used primarily to dis- 
play the structure of the system. The Bond Graph method enables one 
to visualize a system as a set of multi-port devices connected by power 
bonds. As a result of this partitioning of the system the energetic 
interactions across the boundaries are clearly shown. For the case of 
the prototype system under study the basic bond graph representing the 
main power flow is as follows : 



Gas 

Turbine 



N^Altemator 



^ Induction 
Motor 



^ Propeller 



Each component Is a multiport which interacts with the rest of the 
system through the power bonds as shown. Associated with each bond 
are the power variables (T, w or e, i) which determine the power flow 
between the multiports. 

The models required for the multiports shown above can be deter- 
mined in various ways. However, regardless of how the models are de- 
termined, they should satisfy the following conditions: 

1) The time history of the power variables at each port when connected 
to the system should in the simulation match that of the corresponding 
prototype parts over the ranges of variables and frequencies being con- 
sidered. 

2) Each model should be as simple as possible consistent with the ac- 
curacy required for the simulation. Unnecessary complexity will in- 
crease the cost of simulation and will tend to obscure the more impor- 
tant points. More important however is the need to adapt the simulation 
to available computer facilities. Frequently, as in this case, the 
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simulation must be simplified because of limitations of the available 
computing equipment. 

There is no unique method of modeling system components. However, 
any method can be classified as being either analytical or empirical. 
Analytic methods rely on use of the basic engineering disciplines and 
in effect derive a model from first principles. Although in principle 
analytic methods can be used to model any device, these methods become 
extremely difficult to use for very complex components. Empirical 
methods are in effect a transfer function approach to modeling. Given 
a component and information about its performance und a r specified con- 
ditions, it is usually possible to infer a relationship between the in- 
put and output variables. For multiple input-multiple output systems 
this technique may become difficult to use. Also, this method gives no 
information about the internal operation of the component. However, for 
studies where only interactions with the rest of a system are considered 
this method is perfectly adequate. 

In this study the models are determined using both methods. The 
alternator model is based on an analytic derivation whereas the model 
for the gas turbine is derived from performance curves. The models used 
for this simulation are discussed in detail in the following sections. 

The question most frequently asked during this study was "Why simu- 
late a system which has already been built?” The final stage in any sim- 
ulation is validation of the results. Although in most cases the simu- 
lation must be validated without prior knowledge of the prototype, models 
which have been proven with one prototype can be used with great confi- 
dence for a large family of similar systems. In effect, the best 
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documentation of a system model is by comparison with observed perform- 
ance of a prototype. A large part of this study deals with collecting 
sufficient operating data from the prototype system in order to validate 
the models used. 

The prototype was selected because it provided the opportunity to 
model an interconnected system which included mechanical and electrical 
machinery with feed-back control on the energy source. Equally im- 
portant was the ability to instrument the system in order to record 
sufficient operating data to validate the model. The prototype selected 
was the gas turbine emergency generator set and bow propulsion system on 
the USCGC BOUTWELL which is located in Boston, Massachusetts. The pur- 
pose of this study is to simulate the effects of the transients in the 
system caused when the bow thruster motor is started. 
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2 . DESCRIPTION OF PROTOTYPE SYSTEM 

The p rototype system is shown in detail in the bond graph shown 
in Figure 1. The prototype system is defined from an arbitrary choice 
of equipment installed on the ship. The emergency generator can be 
used to supply electrical power to any load on the ship as shown in the 
sketch below. The bow thruster is a trainable propeller which is used 
to provide a directional thrust at the bow when docking the ship. The 
unit is retracted into the ship when not in use. The electrical system 
can be interconnected so that the bow thruster is the only load on the 
emergency generator. The prime mover is a single shaft Solar Saturn 
Gas Turbine. For standard atmospheric conditions the engine is rated 
at 1200 HP at 22,300 RPM. The output is reduced through a two stage 
planetary gear box to 1200 RPM. An engine cutaway view, arrangement of 
reducing gear and engine performance curves are given in Figures 2, 3, 
and 4. 

The engine speed is controlled by a Woodward EG series governor. 
The governor system is shown in the schematic on page 20. The two main 
components are the EGA Control Box and the EGB-2C Actuator. The gover- 
nor can control the engine in either isochronous or droop mode. The 
actuator can also control the speed using the conventional mechanical 
flyweight assembly. For this study the system will be controlled 
electrically and in the isochronous mode. The operation can best be 
described by the transfer functions for the components. The transfer 
functions supplied by Woodward Governor Company are of the following 



form: 
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E ! = 



-1.5 AN, 0.00405 S 



NrE 



s /10 + 1 



REF 



(error signal) 



e 3 _ K (s + a) 
E-]_ + ys + o 



(amplifier section) 



0 _ 33.3 a 2 (s + a 3 ) (EGB-2C Actuator) 

E 3 s( s + a^) 

These signals are shown in the governor schematic. Figure 5. 

The error signal is a combination of a voltage proportional to the 
error in speed and a voltage output from the load sensing network which 
provides lead compensation. The governor response can vary over a wide 
range depending on the values of the following parameters: 

(1) Feedforward gain of the load sensing network (0.00405) 

(2) Gain of the amplifier section (K) 

(3) Damping ratio of the stabilizer section (y) 

(4) Buffer spring stiffness (a-^) 

(5) Needle valve setting ( 013 ) 

The alternator is a General Electric machine with the following 
characteristics : 



Rated 


KVA 


625 


Speed 


1200 RPM 


Rated 


KW 


500 


Winding 


3 Phase 


Rated 


P.F. 


0.8 


Connection 


Y 


Rated 


Volts 


450 


Frequency 


60 Hz 


Rated 


Amps 


801 


Poles 


6 


Voltage regulation 


is provided by 


a static excitation 


system. 



The alternator is connected to the ship’s electrical system in the 
following way. 



-15- 




For this study the system is set up so that the emergency generator is 
supplying only the Bow Thruster load. 

The Bow Thruster Unit is shown in Figures 6 and 7. It is manufactured 
by Schottel of America, Inc. The driving motor is a General Electric three- 
phase wound rotor induction motor with the following ratings: 

Rated Volts 440 Rated Amps 440 

Speed 1180 RPM Torque 1550 lb-ft 

HP 350 Frequency 60 Hz 

The motor performance curves are shown in Figure 11. The motor is started 
with full line voltage applied and additional resistance in the rotor cir- 
cuit. Five propeller speeds are provided by shunting the external re- 
sistors in four steps. The shunting of these resistors is controlled by 
time delay relays. The motor is connected to the propeller through a 
vertical shaft and bevel gears. 

The propeller is enclosed in a nozzle as shown on page 22. The 
diameter is four feet, pitch 28.819 inches and the reduction ratio from 



the motor speed 2.079:1. 



PROTOTYPE SYSTEM BOND GRAPH 
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Figure 2. 




Drawing Courtesy of SOLAR Inc 
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Figure 3 * 
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Figure 5 . 

Governing System Schematic 
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Figure 7. 

Bow Propeller With Nozzle 
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3. INSTRUMENTATION AND RECORDING 

Verification of the system model and simulation by comparison with 
the prototype performance requires records of the system variables during 
operation. The optimum choice for recorded variables are the state varia 
bles of the system model. In this case the choice of observed variables 
was limited to measurements which could be made without major modifica- 
tions to the prototype and which did not require excessive time or cost. 
As a result, most of the observed variables were chosen because the proto 
type was instrumented to monitor these variables. However, these con- 
straints did not prevent obtaining a set of observables sufficient for 
model verification. 

The complete instrumentation setup is shown in the bond graph in 
Figure 8. The main power flow through the system is shown with heavy 
lines. The boxed items show all sensors, signal conditioners and record- 
ers which were added to the prototype instrumentation. Each box is la- 
beled with the variable recorded. All bonds to the boxes are active, i. 
e., zero power flow into the recorders is assumed. Details of each meas- 
urement is given below. 

Results were recorded on a Clevite Mark 220 Brush Recorder. This 
instrument is very versatile and provided satisfactory records of all 
variables measured. Input signals ranged from five millivolts to fifty 
volts and in frequency from D.C. to 60 Hz. For full scale deflections 
the frequency response is flat up to 40 Hz. The frequency response can 
be increased to 100 Hz if the meter movement is limited to twenty percent 
of the full scale. 



ALTERNATOR LOAD The EGA Control Box Load Section senses the 
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alternator terminal voltage and current through potential and current 
transformers. A correction is made for power factor so that the output 
of the load section is a signal with a D.C. level proportional to the KW 
load. At full load the D.C. level is approximately nine volts. However, 
this value will vary depending on the values of the resistors used with 
the current transformers, the levels of voltages and currents delivered 
by the instrument transformers and the gain of the load sensing section. 

The load signal contains harmonics of the line frequency in addi- 
tion to the D.C. level. Since the amplitude of these harmonics was ap- 
proximately twenty percent of the D.C. level an RC low pass filter with a 
time constant of 0.068 seconds was used to filter the harmonics. 

FIELD CURRENT The full load field current is approximately seventy 
amps at seventy volts. A shunt resistor was inserted in the field circuit 
at the brush rigging. The shunt resistor voltage drop is fifty millivolts 
at seventy-five amps. 

ALTERNATOR SPEED The installed tachometers and frequency meter are 
used primarily for steady state observations. Large errors in measuring 
transients can be expected when the time interval of the transient ap- 
proaches the time constants of the meter movements. As a result, the 
meters cannot be used to accurately record transients resulting from step 
changes in load. 

The speed measurement was made using a -Clevite Corporation Frequency 
Converter. The converter is calibrated to give a D.C. signal proportional 
to the change in frequency. A deviation of ± 5 Hz from the base frequency 
gives a ± 2.5 volt output. The rise time of the output signal is 0.060 
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seconds. This response together with the response of the recorder gives 
an accurate measure of the speed transient. 

The signal input to the converter is at 110 volts in the range 55-65 
Hz. The prototype system was set up with the alternator supplying the 
emergency switchboard so that the bow thruster was the only load. In 
order to provide a 110-volt signal at the alternator frequency for input 
to the converter, the 110 -volt lighting circuit at the emergency generator 
was supplied from the emergency board. This total lighting load is less 
than five KW so that the effect on the alternator is negligible. The al- 
ternator load is assumed to be from the bow thruster motor only. 

TERMINAL VOLTAGE The terminal voltage was recorded using a Clevite 
Corporation Volts Converter. The recorder delivers a D.C. voltage propor- 
tional to the peak A.C. voltage. The converter was driven by the 220 volt 
output of the load sensing network potential transformer. 

ARMATURE CURRENT The armature current was recorded from the voltage 
across the input resistors to the load sensing network. The current to 
these resistors is supplied by the current transformers connected across 
the main line. Since this signal was at the 60 Hz line frequency the 
maximum amplitude of the recorded signal was limited to approximately 
twenty percent of the chart deflection in order to have satisfactory fre- 
quency response. 

BOW THRUSTER SPEED The induction motor drives an A.C. tachometer 
generator. The tachometer output was recorded directly on the Brush Re- 
corder. Speed can be determined either from the instantaneous frequency 
or voltage amplitude. For 1200 motor RPM the tachometer frequency is 



i 
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40 Hz. The tachometer signal was recorded full scale on the chart since 
the recorder frequency response is flat to 40 Hz. 



RECORDER 
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Figure 8* 

IIISTRU1-IEI STATION BOND GRAPH 




HJ Load 




-28- 



4.1 GAS TURBINE MODEL The derivation and final structure of the 
gas turbine model depends on the type of study involved. In one type of 
study, knowledge of internal behavior in response to changes in environ- 
ment and load is required. In other studies, only knowledge of the time 
history of the i nput /output variables is required. Obviously for the 
latter case the model can be simpler and will most likely be in the form 
of a transfer function. 

In reference (10) a model is derived from basic principles for a 
single shaft gas turbine where it is shown that the gas turbine can be 
represented by a multiple input-multiple output transfer function of the 



following form. 







The matrix representing the input and output vectors is a relatively com- 
plicated non-linear function of the internal variables and engine charac- 
teristics. For this case the model can be simplified considerably and 
will be used in a linearized form. 

The most valuable information provided by this formulation is an indi 
cation of the functional relationships which exist between the input and 
output variables. Since the mass flow rate (m) is not being considered, 
only the functional relationship for torque is needed. From the transfer 
function it is seen that 



T = *(P lf T ls p 4 , N, mft) 
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For this study, the ambient conditions and can be assumed constant. 
For a constant mass flow rate, fixed geometry and constant ambient condi- 
tions the exit pressure is also relatively constant. Since the mass 
flow rate is primarily a function of the speed, which is relatively con- 
stant, the only input variable which changes significantly is the fuel flow 
rate m^ . As a result the functional relationship can be reduced to: 

t = $> T (mfb) 

This functional relationship can also be derived by considering a 
control surface around the engine. 



Air Flow 
Fuel Flow 



-> 




Air Flow 
Shaft Power 



For this case where the engine can be considered small with respect to its 
ability to store energy as a result of either mass or thermal capacity, 
the changes in stored energy can be neglected. From observations on the 
prototype, it is known that the temperature of the exhaust gas is relative- 
ly constant throughout the operating range. With a constant air flow rate 
and exit velocity the enthalpy of the entering and leaving gases is rela- 
tively constant. As a result, the primary energy quantities entering and 
leaving the control surface are in the fuel flow and shaft work. With all 
of the assumptions above it is seen that output power is proportional to 
input fuel. Since speed is relatively constant it follows that 

t = $ ' ( mf]-, ) = K-rm^ 

This relationship is verified by the engine performance curves given 
in Figure 4. Since the output KW is proportional to torque at constant 
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speed, the torque is proportional to fuel flow rate for constant T^. 

In a simulation, it is necessary to consider the effects of the 
engine dead time or transport delay. This delay is defined as the time 
interval between a change in fuel flow rate and the change in torque out- 
put. From the control surface above, it is seen that this delay depends 
on the ability of the engine to store energy and the speed at which changes 
in internal temperature resulting from fuel rate changes can propogate 
through the engine. As stated above, the storage effects can be neglected. 
An estimate of the delay due to propagation time can be made by consider- 
ing the transit time for sound waves through the engine. This is the speed 
at which the temperature changes resulting from fuel rate changes propogate. 
For the Saturn engine, the transit or dead time is on the order of 0.005 
seconds. The torque relationship is then modified as follows: 



_ _ -Ts* ~ K t m fb 
t - K T e mfb - 1 + Ts 

The delay is approximated by a first order lag for the purposes of analog 
computation. 

All that is necessary to develop the turbine model is to determine the 
torque coefficient from the performance curve. This coefficient is a func- 
tion of ambient temperature with K T decreasing with increasing ambient tem- 
perature. It will be shown below that when the system block diagram is 
non-dimens ionali zed for the analog computation, the torque constant can be 
lumped with the overall loop gain. 

4.2 GOVERNOR SECTION The governor schematic is shown in Figure 5. 



The engine block diagram can be formulated directly from the schematic us- 
ing the governor transfer functions supplied by Woodward Governor Company. 
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Figure 12. 



From the diagram it is seen that if the delay is neglected 
T e = K'Aa = K T ^ f 0 = K T »K f 0 

where it is assumed that the fuel flow is proportional to the angle of the 
fuel control valve. Although in most cases this relationship is approxi- 
mately linear, slight non-linearities will not affect the overall system 
performance. Although the fuel control valve is actually a dynamic opera- 
tor, in this case the valve will be considered to be a static operator. 

For low frequency changes in fuel angle this is a good approximation. 

With K t = K T ! Kf ; t= K t 0 Taking t = t Ref when 0 = (0 max - O m in)> 

where (0 max - 0 m i n ) t ^ e roaximum travel of the fuel control actuator, 

the torque constant K is TD rr 

1 v - REF 

K x " 0 _ q . 

u max u min 
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In reference (10) it is shown that if the loop gain 1<2 is not exces- „ 
sive the gain and phase margins at the crossover point of the open loop 
transfer function (approximately 6 rad/sec) will be adequate and the action 

of the stabilizer section of the control box f ) can be neglected. 

This simplification neglects the response of the governor system to fre- 
quencies above approximately 50 rad/sec. With the stabilizer section ne- 
glected and using the nominal values for the governor settings, the system 
block diagram is as shown in the section on Analog Scaling page 58. 

4.3 ALTERNATOR MODEL The study of alternators typically involves 
deriving an equivalent circuit and phasor diagram of the electrical varia- 
bles. However, the phasor diagram does not explicitly show the mechanical- 
electrical energy conversion process which is a major part of the prototype 
system’s power interactions. The model used for the simulation must in- 
clude the torque-current and speed-voltage relationships involved in the 
energy conversion process. Such a model can be derived from a study of 
the energy interactions of idealized magnetic circuits. The model proposed 
by Paynter, reference (10), was formulated this way and as a result shows 
the structure of alternator as a gyrator-type converter as well as the 
behavior of the electrical variables. 

The major assumptions required to simplify the analytical derivation 
of the alternator model are that the magnetic circuits are not saturated 
and that all fluxes are linear functions of the currents. The lineariza- 
tion reduces the problem to a study of linear-coupled circuits and allows 
superposition of the magnetic fields. It is also assumed that the air 
gap is uniform. However, in a study of terminal voltage, current, torque. 
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and speed relations under load conditions, salient pole effects are rel- 
atively minor and can be neglected. One of the main effects of saliency 
is to produce a reluctance torque, and, except at no load, this torque 
is negligible compared to the total torque. The resulting model can then 
be used for machines which are operated under quasi-linear conditions and 
can be used for modeling the prototype alternator. 

The basic variables and relations required for the derivation are as 
follows : 

A 5 flux linkage = N4> 

N = effective turns ratio of winding 

4> a equivalent flux linking all terms = PF 

P a permeance 

F a magnetomotive force = Ni 

F i 

W m * a co-energy in magnetic field = / 4>(F)df = Jx(i)di 

o o 

The inherent behavior of polyphase machines is the same regardless of 
the number of phases. As a result the torque is determined by the current 
flowing in each phase and the number of phases. In the balanced polyphase 
machine each phase torque is a pulsating torque at twice the line fre- 
quency with a D.C. level which is exactly the total torque divided by the 
number of phases. When all phase torques are summed the pulsating torques 
cancel so that the total torque is constant in the steady state. 

A reasonably accurate model of the alternator (and motor) can be made 
using only a single phase if the torque pulses can be neglected. Alternators 
with large inertias behave dynamically as low pass filters as seen from the 
simplified bond graph and system equations for a rotating mass. 
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Where — = — - Where I = 

T Is + R g 

For the prototype system I is on the order of 2000 lb-ft 2 and 2co 

the frequency of the torque pulsations is 754 radians per second. Since 

R is relatively small; 



0) 

T 



1 



~ 2000 
32 



x 754 



<< 1 



where oo is the fluctuation in speed 
caused by the pulsating torque. 



The inertia of the induction motor and load is the same order of mag- 
nitude so that a single phase model can be used for both with torque pulsa- 
tions being filtered by the relatively large inertias. 

The alternator model used below is taken from one phase of a two- 
phase alternator model. Although the model is derived for a two-phase 
machine, one phase can be used to represent one phase of a three-phase 
machine since the torque pulsations are both double -frequency pulsations 
and are filtered by the inertias. All that is required to use this model 
for a three-phase machine is to properly scale the D.C. level of the 
torque . 

The two-phase representation shown below is very similar to the Park’s 
dqo transformation used in studying three-phase machines. In this case the 
variables are referred to the two-phase stator whereas in the dqo trans- 
formation the variables are viewed with respect to the rotor. 
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Total flux linking coil a 
*a = N a P[N a i a + N f i f C] 

Total flux linking coil b 

*b = N a P[N a i b + N f i f S] 

Total flux linking coil f 

X f = N f P[N a i a C + N a i b S + N f i f ] 

Where N a = N b M 5 N a NfP 

C = cos 0 N = N a 2 

S = sin 0 



( 1 ) 



( 2 ) 



(3) 



and the permeance is independent of 0 in accordance with the assumption 
of a smooth air gap. 

Now W m * = /[A a di a + A b di b + A f di f ] 

W m * = J[PNi a di a + CMifdi a + PNi b di b + SMi b di b + CHi a dif + SMi b di.f 
+PN f 2 ifdi f ] = ^PN[i a 2 + i b 2 ] +^PN f 2 i f 2 + Mi f [Ci a + Si b ] 



( 4 ) 
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For any electro-mechanical converter 

= T S ° T = ^f E“Si a + Cibl 

The exact relationship between t, w, e and i is found from a con- 
sideration of power balance. 

Assuming a lossless system: 
p mech = P AG (developed at air gap) 

P elect output ~ P AG p stored magnetic field) 

The alternator is in effect separated into two processes, energy conver- 
sion and energy storage. Equating power at the air gap yields the fol- 
lowing identity: 

™ = P AG 

Where for a two-phase machine 
to) = e a i a + e b i b 
MifW C-Si a + Ci b ] = e a i a + e b i b 
or e a = [-SMif] w = Ra tu 
eb - [CMif ] a) = 00 

Where R a = - SMif and Rb “ CMif 

in matrix notation 



e a 




1 

O 

o 

_°L_. 




i 3 


e b 


= 


p? 

O 

o 


X 


L i> 


T 




R a R b 0 







If x is separated into its two components the gyrator-like action of the 
energy conversion is evident. 

T a = Va e a = R a“> 

Tb = R bib e b = R b w 



- 37 - 



The magnetic energy can be modeled as an I field. The three currents 
i a , ifc , and if are the sources of the stored magnetic energy. Since 
if << i a and 

if << ib 

most of the energy is stored as a result of the armature currents so that 
changes in armature currents will account for most of the change in stored 
energy. This implies a model of the form 

I 

e a v t 

ia " 1 "ia 

which is in fact the equivalent circuit of an alternator with negligible 
armature resistance. 



O 'boo&bboo&zr* 



The bond graph for the complete conversion process is: 

I, 




N 1 



H 1 h 



V4 



X b 



VFO = variable 

frequency 

oscillator 






Z 
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where E represents the prime mover. The causality of E is determined by 
the structure of the complete system. The Analog simulation of this 
model will represent the exact alternator behavior (within the limits of 
original assumptions). However , the model requires six multipliers. The 
need to minimize Analog hardware required for simulation makes it neces- 
sary to use the single-phase model. 

For this case 

T a = Va 

= -MIf sin ( w t) i a 

for i a = I a sin(a 3 t -<j>) where 0 = cat 

<j> 5 phase angle between i a and e a 

t = - MIfI a sin (cot) sin (cot -$) 

= MIfI a sin (cot) sin (tot - <f>) 

= ^MIfI a [cos<j> - cos ( 2 cot - <f>)] 

Whereas T a + T b = ^f^a cos ^ » 

T a alone gives a D.C. value equal to \ of the actual steady torque 
plus a double frequency pulsating value. As shown above the alternator 
rotating mass has a large amount of inertia and acts as a filter for the 
high frequency torque pulsations. As a result, the one phase approxima- 
tion can be used without any appreciable loss in accuracy. 

The total torque input to the prime mover can then be obtained merely 
by doubling the torque calculated above. 

For the simplified model the bond graph becomes: 
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/ 








In this model the inertia of the rotating mass has been included and is 
shown connected to the -1 junction of the alternator. Now the causality 
required for the prime mover can be determined. 

The I element imposes the velocity at the 1 junction so that the fol- 
lowing causality must be imposed. 




MGY 



This requires that the prime mover be modeled as some type of torque 
source. This is as expected because if a velocity source were assumed, 
step changes in the input would impose a step change in the velocity of 
the rotating mass and to do so would require infinite power. As a result 
the model of the gas turbine must be of a form such that velocity is the 



-40- 



input signal and torque is the output signal. 

From the simplified model bond graph the state variable representa- 
tion can be formulated from the following equations. 



d<% 

dt 

di - 
dt 




s 



“ T l } 
- v) 



e = 

t = Ri 

R = MIf cos a)t 
v = Zi 



This model involves five parameters. 

I E inertia of rotating masses 
eq 

= synchronous reactance of alternator 
M e this parameter is defined above as N a N^P. However, it is 
closely approximated by the slope of the air gap voltage 
line of the alternator at synchronous speed, 
e = [MIf(cos wt)] w 




constant 



If E field current, for the initial simulation will be assumed 
constant . 

Z E connected load. 

The Analog setup sheet is shown on page 57 where the load Z is re- 
placed by the induction motor model. 
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4.4 INDUCTION MOTOR A commonly used induction motor equivalent 
circuit is the Steinmetz equivalent circuit shown below. 



-WA/W TrsThsy- 



Ri 



X, 



Xo 



^m 



:Ws 



The derivation of this equivalent circuit is given in standard Electrical 
Engineering texts. This model can be simplified considerably without 
serious error by neglecting the magnetizing reactance X m . The circuit 
can then be simplified to the following. 



o r 6docd3r- 

^eq 



X eq = X x + X 2 



>V S 



R eq = R 1 * 51. and R. 



Except at standstill Req ^ R 2 /s. In the case of the wound rotor induc- 
tion motor large resistance is added to the rotor circuit for starting so 

that through the entire operating range the equivalent resistance is 
R 2 ? 

R eq % ~ wherp V = R 2 + R added • 

For the case of the prototype motor, the added resistance is approxi- 
mately thirty times as large as R 2 and is shunted out in four steps so 
that Rq can be neglected at all speeds. 

By applying a power balance across the air gap and neglecting losses 
an expression for the developed torque can be derived. 



t = ' qi? 
m "s s 



where ( w s ) is the synchronous mechanical speed and (q) is the number of 
phases . 
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4.5 PROPELLER The torque relationship for any non-cavitating 
propeller can be represented by an equation of the form 
x = p D 3 n 2 C T (g) 

where the water density (p) and diameter (D) are constants. The torque 
coefficient C X (B) is a function of the velocity angle (3) where 3 is de- 
fined as 3 = tan“^(V^/V w ) where Vq is the axial velocity of flow through 
the propeller and 

Voj is the tangential velocity of the propeller at the maximum diameter. 
The angle 3 is identical to the propeller advance ratio (J) used in propel- 
ler design by naval architects where J = Vo/irnD. With 3 constant the 
propeller can be modeled by an expression of the form 
x -= K p n 2 

From the picture of the propeller shown in Figure 7, it is seen that 
the propeller operates in a nozzle with an axial length approximately one- 
half diameter. With this nozzle the flow rate through the propeller is 
of the approximate form <*Q/D 

V * n 

As a result the velocity angle 3 becomes 

3 = tan -3 -(-2— ) - const 
TinD 

The torque expression given above then reduces to the form commonly 
used to represent propeller loads over a wide range of RPM where torque 
is taken as directly proportional to RPM squared. 



The prototype system was operated with the ship tied to the pier 
when the data was taken. Although there is some slight motion of the 
ship during the operation of the bow thruster it was assumed that the 



-43- 



propeller is stationary and that the axial velocity of water through the 
propeller is due to water motion only. 

4.6 BEARING RESISTANCES All bearing resistances will be neglected 
in the simulation. In the case of the alternator where the speed devia- 
tions are slight the bearing drag is constant and will have no significant 
influence. For the induction motor the resistance will vary, but the 
resistance effects are assumed to be small compared to the load and driving 
torques. The resistances can therefore be neglected with no serious error. 
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5.1 PARAMETER ESTIMATES The system model to be used for this 
simulation requires values for the following parameters: 

1 ^ - engine inertia 

Ig - accessory drive inertia 

1 3 - speed reducer inertia 

1 4 - alternator inertia 

X sa - synchronous reactance of the alternator 

X sm “ synchronous reactance of the induction motor 

R a - rotor resistance of motor 

R a ! - added series rotor resistances 

I - motor inertia 



K t - engine torque/fuel flow ratio 
governor parameters 

Where possible these parameters were obtained from information sup- 
plied by the equipment manufacturers. However, in some cases it was 
necessary to estimate parameters. 

I]_ The engine inertia is divided into three parts, the compressor, 
turbine and accessory drive package. Estimates for compressor and turbine 
inertias were obtained from Solar. 



m 



Compressor 



1.750 in-lb-sec^ 



Turbine Stage (1) 



.532 



Turbine Stage (2) 



.564 



Turbine Stage (3) 



.916 



3.762 in-lb-sec 2 at 22,300 RPM 
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The accessory drive is approximated as follows. The idler gear 
and output gear are the largest rotating parts at the intermediate speed. 
These are approximated as steel discs one-half inch thick and eight inches 
in diameter. For a single disc, 

I •= kmr 2 



m = 



TTD ^TP 
49 



where for steel P = 487 lb/ft 



m = tt ( 6 4 ) 487 = 0.0184 lb -sec 2 /in 

4(2) 1728-384 

I = ^mr 2 = 0.147 in-lb-sec 2 

For two discs and considering the remaining parts as having the same inertia 
as a single disc, the total inertia of the accessory drive is approximately 
0.5 in-lb-sec^ at 6010 RPM. 

*3 The speed reducer inertia is estimated by assuming that all 



inertias are negligible except that of the first-stage ring gear. Since 
all inertias will be referred to the output speed, both speed ratios and 
WR^ values must be considered. The first-stage ring gear is approximated 
by a one-quarter inch steel disc eleven inches in diameter plus a ring two 
inches wide by one-quarter inch thick and eleven inches in diameter. The 
resulting inertia is approximately 1.0 in-lb-sec^. Although the inertias 
of the second-stage parts are of the same order of magnitude as the first 
stage ring gear, they become insignificant compared to the first stage 
parts when referred to the output speed. 

-^4 The alternator rotor inertia was given on the manufacturer’s 



I = 29.6 lb-ft-sec 2 



drawings . 

WR 2 = 951 lb-ft 2 



or 
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All inertias can now be referred to the output speed of 1200 RPM. 



X lEQ = I 1 



fN 



turb 



N 



3.76 

12 



ait 

( 22 ,30C ~ 
1200 



■EQ 



= I r 



rN 



acc drive 



N 



alt 



= 108 lb-ft-sec^ 
2 



0.50 r6010)2 



12 



X 3EQ “ X 3 



I | - 1 lb-ft-sec 

^ 1 2 0 0 j 

M, .2 

1 1st stage 



'•alt 



1.0 r6010l 2 



12 <1200 



2 lb-ft-sec^ 
2 



I 4e q = I 4 = 29.6 lb-ft-sec 

It is seen that the inertias of the engine and the alternator are the 

2 

dominant terms. The total inertia is taken -as 140 lb-ft-sec . 

^sa The alternator synchronous reactance is given cn General 
Electric drawings as 1.19 ohms. For a line frequency of 60 Hz; 

y 

L = - s — = 0 . 00316h 
0) 

^sm, Fa The motor reactance can be estimated from the equivalent 
circuit model used for the motor and the motor performance curves. 

N s -N 



X, 



‘sm 



r 9 

R = — 



•s = 



Nc 



From performance curves at the motor rating point: 
N - 1180 RPM 

I = 440 amps s = 1.66% 



E = 440 volts 
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For the equivalent circuit used, the magnitudes of the phasor quantities 
are related by: 

E . 2 

I = - where Z = (R + X 2 ) 2 
Z s x 

X s = 0.71 ohm 
E 

Z - ~ - 1 ohm and L - 0.00188h 

The rotor resistance is given in the General Electric drawings as 0.0117 ohms 
^•m No figures were available in the manufacturer's drawings for the 
motor inertia. The inertia was estimated from a drawing of the rotor where 
the mass was approximated by a steel right circular cylinder fifteen inches 
long by twenty inches in diameter. The inertia estimate is 14.1 lb-ft-sec^. 
This is approximately fifty percent the inertia of the alternator and is a 
reasonable estimate. 

The propeller is approximated by a steel disc four feet in diameter 
and one inch thick. This gives a value of 30 lb-ft-sec 2 . When this iner- 
tia is referred to the motor speed through a speed ratio of approximately 
2:1 the equivalent inertia of the propeller becomes 7.5 lb-ft-sec . 

The ratio of steady state engine torque to fuel flow is obtained 
from the engine performance curves. For constant speed and constant ambient 
conditions the torque is linear with fuel flow throughout the operating 
range. In this study where fuel input is not observed it is not necessary 
to isolate this constant from the overall gain in the engine control loop. 

In scaling of the simulation it is shown that where torque and fuel actu- 
ator angle are considered in non-dimensional form this constant is reduced 
to unity. 
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No information on the propeller was available from the manu- 
facturer. However, for the assumed propeller model which is 
tp - V 2 

a torque coefficient can be fit from the observed propeller performance. 
Using the maximum observed propeller speed for fitting the coefficient 
gives the following: N = 1125 RPM 

Neglecting frictional losses the motor torque at this speed is 1370 lh-ft 
so that 



"p-S* 

= 0.00108 lb-ft/RPM 2 

GOVERNOR SETTINGS There are five parameters connected with the 
engine governing system. These parameters can vary over a wide range 
depending on the governor settings. In order to determine the exact 
settings it would be necessary to measure the open loop frequency re- 
sponse of the control system. The ranges of the governor settings as 
well as nominal values were supplied by the Woodward Governor Company. 
The nominal values will be used in this study. The feedforward gain and 
loop gain will be varied and the effect on system performance observed. 
The values for the other parameters will be taken as follows : 

Parameter Range Nominal Value 

23 to 184 



“1 

c»2 



28 

1.6 



“3 

K 



0.08 to 160 
vary in simulation 



4 

30,000 



a 



10 
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Parameter Range 

Y 
a 

Feed Forward Gain vary in simulation 

R a The added rotor resistances are given in the 
drawing. The resistances are shunted out in four steps 
Step 1 .240 ohms /phase 

Step 2 .088 ohms /phase 

Step 3 .072 ohms /phase 

.063 ohms/phase 



Nominal Value 
46 

2210 

0.0045 

manufacturer ' s 



Step 4 
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5.2 SIMULATION-ANALOG SET UP The models derived in previous sec- 
tions can be set up for Analog computation in a straightforward manner. 
The set up for the complete system is shown in Figure 9 on page 57. 

GAS TURBINE The Analog set up for the gas turbine and speed control 
is implemented directly from the block diagram of the engine on page 31. 

(1) The feedforward of load disturbance is a high pass filter of the 

form: 



Y = KS 
X TS + 1 

This transfer function is realized by a circuit of the form below 
which is taken from reference (3). 




where 



A = 1/T 
B = k/T 



The operator is realized using A40 and A31 in the simulation. 

(2) The engine speed is obtained at A60 by integrating the difference 
of the driving and load torques. 



A 





x a ) dt 



0 l( ) - ce (t Q ) 



Since co (t Q ) can be taken as the base speed of the alternator the 
output of A60 when properly scaled gives the error in engine speed. This 
voltage is summed with the output of the feedforward network to give the 
error signal developed in the EGA control box. 

(3) The error signal is generated by A61. 



Kl s 



E -i — x p ~ Ka 

1 T s + 1 a 2 <u 



Ace 



where 



ref 



T = 

a 



alt 



t ref 
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(4) Since the operator representing the stabilizer section of the 
control box has been neglected the transfer function between the error 
signal and the output of the fuel actuator • angle 0 is 

_ K 3 (T t s + 1)(T 2 s t 1) 

E 1 S (T 3 S + 1) 

This transfer function can be realized by taking two operators cas- 
caded . 

Y Tis +1 

For — = — 

X s 

Y = T X X t X 
or Y = T L X t /X 

This operator is realized using A95 , A39 , and A01 




Y (T 9 s +1) 

The operator for — = - — * — is also from reference (3). 



X (T 3 s t 1) 




A = K t 2 /T 3 
B = (T 2 -T 3 )/T 2 T 3 
C = i/t 3 



The output of A90 represents the angle of the fuel control actuator. In 
the derivation of the turbine model it was shown that except for a trans- 
port delay the torque is proportional to fuel input. The fuel control 
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angle is assumed proportional to fuel input. The engine dead time can 
be neglected so that except for a gain factor the output of A90 represents 
the gas turbine torque. 

With the correct setting on P60 the output of A91 represents engine 
torque . 

The gas turbine model is now complete except for the load torque. The 
alternator torque is generated bv the multiplier A38. This torque is of 
the form 

T = T q (1 + COS 2 Wt) 

It is necessary to filter the double frequency component for two reasons. 

(1) The pulsations are- a form of noise which cause slight variations 
in engine speed. When compared to the base engine speed of 1200 RPM, 
these fluctuations in speed are insignificant. However, with a fine speed 
control on the engine where the variations in speed due to the real load 
changes are of the order of 10 RPM, the noise from the load torque signal 
becomes significant. This point is most evident when scaling the integrator 
which generates the speed signal. If the integrator is scaled so that the 
output is N/N ref where N re f is 1200 RPM the integrator gain which represents 
the reciprocal of inertia is small and corresponds to a large inertia. In 
the simulation it is necessary to generate N/N T where N’ is of the order of 
the maximum deviation in engine speed due to changes in load torque. As a 
result, the integrator gain is increased to the point where the double 
frequency torque pulsations are passed along with the D.C. torque level. 

(2) It is also necessary to filter the load torque because it is in- 
put to the load sensing network which has a derivative action. 
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A low pass filter was chosen over a band stop filter for removing 
the torque pulsations because the low pass filter requires less Analog 
hardware. An attenuation of 40 db at the noise frequency will provide 
adequate filtering. At this level the noise is barely noticeable on 
the oscilloscope. A first order filter tuned for an attenuation of 40 
db at the noise frequency of 120 Hz will have its break point at 1.2 Hz 
and is considered to be too low for the expected load transients. As a 
result a second order filter was used with the break point at 12 Hz, and 
40 db attenuation at 120 Hz. The transfer function of the second order 
filter is of the following form: 



G( jw) 



1 



1 + 2 C 


isi 




'ju| ‘ 




^ n . 


+ 


,'wn J 



With the break point at 12 Hz and with C set at 0.55 in order to minimize 
overshoot near the break frequency, the transfer function becomes: 



G(s) = 



1.1 _2 
1 + 7~5 S + -4 

/-b 75 2 



where s = jw 



Performing the operations indicated gives 
Y = 5650X - 92.5 Y -5650Y. 

Scaling this equation to a form suitable for Analog computation gives 

- 56. 5Y. 



— 




’ 


Y 




Y 




= 56. 5X - 82.5 




100 




100 



This filter is shown on the Analog set up sheet using A70, A74 and A100. 

ALTERNATOR, INDUCTION MOTOR, PROPELLER The equations representing 
these components were derived previously. 
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(1) ^ = i <t e 

dt I 2 e 






(2) R = MIf cos wt 

(3) x a = R. 

d n * i 

(4) 



dt (L a + L m ) 

(5) e = Rw 

(6) v 

(7) x 



(e - v) 



•H 

II 

CM 

•H 

II 


r 2 0) 0 


i R 2 


Nn 


S 


0)0 -w 




1 s; 
o 

1 

\ ^ 



m 



i R 2 



= l'V 



(8) T p — Kr 



(9) 77 = 



dN 

dt 



30 



(Ip + 



(x m - y 



Equation (1) has been described above with A60 generating the 
alternator speed. 

Equation (2) determines the gyrator modulus of the alternator and 
is a function of the field current and frequency. In the prototype the 
speed variations were insignificant compared to the base speed so that 
the frequency can be considered constant. Although the field current 
did vary slightly a constant field current will be used. The gyrator 
modulus is then of the form 

R = K cos (tut) 

The 60 Hz sinusoidal waveform was generated using the Wavetek signal 
generator. 

Equation (4) is solved using A00 where e is the air gap voltage of 
the alternator and v is the back emf of the induction motor. 
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Equation (5) is the air gap voltage of the alternator. With the 
assumption of constant frequency and field current, e is a sinusoidal 
voltage of constant frequency and amplitude. 

Equation (6) generates the back emf of the induction motor and is 
realized at multiplier A03 . 

Equation (7) generates the torque output of the induction motor. 

This torque also contains a double frequency component. However, since 
the induction motor speed is scaled to the 1200 RPM synchronous speed 
it is not necessary to filter this torque. 

Equation (8) is realized at multiplier A63 where the propeller 
torque coefficient is set on P67. 

The switching of the induction motor rotor resistances is set up 
in the network which is input to the junction inverter A04. The values 
of the switched resistors are set on the potentiometers P95, POO, P32, 
and P65 with the residual resistance on P03. The switched resistors are 
input to the amplifier junction of A04 through the D/A switches shown. 

The switching of the D/A switches is controlled by the logic outputs of 
the comparators . 

In the prototype system the switching is controlled by time delay 
relays. This switching action is simulated in the following way. AOS 
generates a ramp timing signal. This signal is input to comparators C04, 
C34, C64, and C94. These comparators have a negative bias corresponding 
to the time at which switching is to occur. The inverted logic output 
of each comparator is input to the switch control of the respective D/A 
switch. Initially the output of A05 is zero. All comparators are in the 
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low position. As a result, the inverted output to the D/A switch is high 
and all switches are conducting. The input to A04 is 
1 ( + R]_ + R 2 + R 3 + Rq. ) . 

As the output of A05 increases to a level above the bias of each comparator, 
switching occurs when the comparator logic output goes high and the in- 
verted output goes low so that the respective D/A switches are turned off. 
After all switching occurs, only the resistance due to P03 remains in the 



circuit . 
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5.3 SCALING ANALOG SIMULATION 

GAS TURBINE When nominal values are substituted for the governor 
settings the block diagram for the prime mover and governor is as shown 
below. 




The machine time constant (T ) is defined from the non-dimensional equa- 

m 

tion for the alternator speed. 



T e - T a - 




d 


r u ' 


T 0 


l T oJ 


dt 


o 

3 



IU3n 

where T = — r— 
m T n 



The base speed o)q is 125.6 rad/sec and the inertia I is the combined iner- 
tia of the engine 5 speed reducer and alternator referred to o)q . The selec- 
tion of the base torque is somewhat arbitrary. Based on the engine power 
rating of 1200 HP with output at 1200 RPM the base torque is: 



1200 x 550 _ - 

T 0 ~ 125.6 " 5250 lb ft 



The maximum expected steady load torque is approximately 1500 lb-ft. The 
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maximum engine torque is limited by the acceleration limiting response 
of the fuel control valve. For a step load change of 1500 lb-ft the 
maximum engine torque developed is assumed to be 3000 lb-ft. With this 
base torque the engine time constant is 5.86 sec. The corresponding gain 
on the integrator which generates Aw is 0.170. The expected maximum 
deviation in speed is taken to be 7.5 rad/sec. With this base speed the 
gain on the speed integrator is: 



Aw 

7.5 





a_ 



where T_ T 
m 



140 x 7.50 
3000 



0 . 352 sec 



The gain on the speed integrator must be increased by a factor of 16.7 
for this new value of of T m T . This is accomplished by modifying the sys- 
tem transfer function as shown below. 
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The operator , is generated by A60. The scaled equation is* 
■‘m s 



f A “"| _ 10(2840) 


~ T e - T a] 


L 7 * 5 J 


3000 



with the pot settings on P38 and P97. 

The scaled feedforward transfer function with T = 10 and K = 0.1 is 

generated at A40 and A31. The gain of the feedforward is then set on Q07 

The operator — + - is generated at A01. For T. = 4 the setting 
s I 

on P08 is 2500. 

^ K(s/10 + 1) . 

The operator 1S generated at A90. In order to reduce 

the number of pots required K is taken as IOOT3/T2. With T2 = 1/10 and 
T3 = 1/28 the pot settings and gains are: 

P62 = 100(2800) 

P37 = 100(1800) 

The inputs to A65 and A91 from A01 are then at a gain of 100. Since the 

gain of 100 is difficult to obtain at the input to summer A91 the output 

of A91 is scaled down by a factor of 10. The inputs and output of A91 

are then as shown on the Analog diagram, Figure 9. The output of A91 cor 
10 T3 

responds to — 0. 

t 2 

0 

The loop gain is introduced at A90. With K = q _ -q \ — and 

max “ u min 

®max “ ®min taken as 30° the required gain at P60 is determined as shown 

* The following convention is used to represent all scaled equations: 
normalized variables in machine units are in brackets, potentiometer co- 
efficients are in parentheses and amplifier gains before the parenthesis. 
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below. 



92.5 K-rO 



92.5Q 

®max “ ®min 



The output of A91 is 10 -2- 0 



10 T 3 (gain) 0 = 92.5 Q 



®max “®min 



(gain) = 



92.5 T 2 



= .8650 



T 3^®max " ®min) 

The required setting on P60 is then 1(8650). 

LOAD COMPONENTS The induction motor speed is generated at A35. 
The unsealed equation is: ^ 



N = 



Is ^ Tjn T P^ 



The maximum motor speed is 1200 RPM which is the synchronous speed. The 
maximum steady motor torque is approximately 1500 lb-ft. Since the simu- 
lated motor torque has the double frequency component the maximum torque 
is twice the instantaneous D.C. level. As a result the maximum torque 
is taken as 4000 lb-ft. With the inertia (I) taken as 21.6 lb-ft-sec^ 
the scaled equation becomes: 



pi 


10(1475) | 


r T m - T p 


|_120oJ 


1 ' ■ 1 


I 4000 



Set on P35 

The propeller torque is generated at A63. The unsealed equation is: 

T p = KpN 2 

_3 

With Kp = 1.08 x 10 and the torque and speed base values as above the 



scaled equation becomes: 
- 1(3880) 



pi 

[4000 I 



f 1 

[l . 44 x 106J 



Set on P67 
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The voltages and current in the simulation must be scaled based on 
peak values rather than the RMS values. For maximum RMS values of 440 
volts and I/phase of 1000 amps the base values are taken as: 

T base = 1414 arr) P s 

^base = 750 volts 

Multiplier A33 generates the unsealed product of e and i. This is 
scaled to the induction motor torque as follows. For the induction motor 
for three phases : 

P(t) = 3 V max sin(wt) I max sin(a)t + $) 
or 

= 3 __ v max^max ^ + cos2qjt:) 

Tm 2^ s 

where the phase angle can be neglected in scaling the equation and (.7376) 
is a conversion from newton -meters to lb -ft. oj s is the synchronous me- 
chanical speed which for this case is 125.6 rad/sec. With the base vari- 
ables assumed above the scaled equation becomes: 



[woo] = 10(2: 340 ) [vro] 



Set on P36 

The output of A00 must be scaled to I max = 1414. The unsealed equa- 



tion for A00 is : 



1 = — 
s 



1 (e -v ) 



Using the base values for current and voltage and the equivalent inductance 
(L) for a single phase the scaled equation is: 

10 3 (8460) 









The motor voltage is generated at A03. The unsealed equation is: 



RaIN s _ RI 
N s -N = 1-N/N S 



v = 
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The denominator is generated at A06 and is already in a non-dimensional 
form. Scaling the equation for current and voltage gives: 



r v n 


1.89 R,, 1 


r 


i — 

<i 

cn 


1 - N/N s 


t 141*+ j 



The potentiometers in the switching network are then set to 1.89 R a and 
R a r where R a is the value for the rotor resistance and R a T the value of 
the added rotor resistance. 

The alternator torque is generated at multiplier A3 8. As shown 
above this output can be scaled to a base torque of 4000 lb- ft using the 
same setting as on P36 which is 10(2340). However, the alternator torque 
is filtered before being input to the gas turbine section so that it is 
not -necessary to scale this signal for the double frequency pulsation. With 
the turbine base torque taken as 3000 lb-ft the gain on the output of A38 
is 10(3100). Since the second order filter is set to a gain of (l) this 
gain is set on P07. 

The scaled variables for all amplifiers and settings for all poten- 
tiometers are given on the assignment sheets given in Tables 1 and 2. 
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Table 1. 

AMPLIFIER ASSIGNMENTS 



Amplifier 


Scaled Value 


Comment 


00 


i/1414 




01 




Governor Control Box 


02 


e/750 


Alternator Modulus 


03 


v/750 


I.M. Back emf 


04 


1 . 89R a i/1414 


For I.M. Back emf 


05 




Time Signal for Resistor 
Network 


06 


(Ng~N)/N 


I.M. Slip 


07 


x a /3000 




09 


v/750 


Invert I.M. Back emf 


10 


(v-e )/750 




14 


-(N s -N)/N 


Invert A06 


30 




LPF for Viewing Torques 


31 




FFWD of Governor 


32 


-i/1414 




33 


vi/750 x 1414 


I.M. Torque 


34 


-N 2 /1.44 x 10 6 




35 


N/1200 


I.M. Speed 


36 


(T p -t m )/4000 




38 




Load Torque on G.T. 


39 




Governor Control Box 



40 



Governor FFWD Net 



ipli 

44 

60 

61 

63 

65 

70 

74 

90 

91 

95 

100 



Scaled Value 



Aco/7.5 
N 2 /l. 44 x 10 6 



T e /3000 
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Comment 
Invert A02 

Change in Alternator Speed 
Governor Error Signal 
For Propeller Torque 
Governor Control Box 
Butterworth Filter 
Butterworth Filter 
G . T . Torque 
Fuel Control Angle 
Governor Control Box 
Butterworth Filter 



00 

02 

03 

04 

05 

07 

08 

10 

32 

34 

35 

36 

37 

38 

60 

61 

62 

64 

65 

67 

70 

90 

95 

97 

100 
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Table 2. 



POTENTIOMETER ASSIGNMENTS 



PARAMETER 

Aw/uq Gain 
FFWD gain 
R a (switched) 

Alt . modulus 
^a 

Switching time 
Timing signal 
Scale x a 

Cont. box setting 
Motor inductance 
R a (switched) 

Switching time 
Inertia of I.M. 

Scale motor torque 
Govt. Cont. box 
Alt. inertia 
Control loop gain 
Butterworth filter 
Govt . Cont . box 
Switching time 
R a (switched) 

Scale propeller torque 
Butterworth filter 
Switching time 
R a (switched) 

Alt. inertia 
Butterworth filter 



SETTING W/GAIN 

1000 ( 10 ) 

2000 (vary) 

6000 

Set for A07 output =.8500 

0597 

1178 

0500 

3103 (10) 

2287 

8277 (1000) 

2000 

3291 

2000 ( 10 ) 

2341 (10) 

1800 (10) 

2841 (10) 

2177 (10) (vary) 

5649 (100) 

2661 (100) 

4895 
2005 
3880 
8243 
6 819 

2001 

2841 (10) 

5648 (100) 
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6. RESULTS 

The results of monitoring the prototype performance and the results 
of the simulation are shown on pages 69 through 75 . A chart speed of 
twenty-five millimeters/second was used for all recordings. Strips A-L 
correspond to a time duration of approximately fourteen seconds while 
channels M-BB span approximately the first eight seconds of the simulation. 
The records are presented in the following order. 

A. Observed Alternator Speed Deviation 

B. Simulated Alternator Speed Deviation 

C. Observed K.W. Load (through 0.068 sec. LPF) 

D. Simulated Alternator Torque (through 0.068 sec. LPF) 

E. Observed Line Current 

F. Simulated Line Current 

G. Observed Induction Motor Speed (A.C. tachometer) 

H. Simulated Induction Motor Speed 

I. Simulated Motor Torque (through 0.068 sec. LPF) 

J. Simulated Alternator Torque (through Butterworth filter) 

K. Simulated Engine Torque - No Dead Time, No FFWD 

L. Simulated Engine Torque - No Dead Time, With FFWD 

M. Alternator Speed with FFWD Gain (Q07) at .1000 

N. Alternator Speed with FFWD Gain at .2000 

O. Alternator Speed with FFWD Gain at .3000 

P. Alternator Speed with FFWD Gain at .5000 
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Q. Simulated Output of Governor Load Sensing Section 

R. Simulated Governor Error Signal 

S. Observed Alternator Terminal Voltage 

U. Simulated Alternator Speed with 0.1 Sec. Dead Time 

V. Simulated Engine Torque with 0.1 Sec. Dead Time 

W. Simulated Fuel Control Angle with 0.1 Sec. Dead Time 

X. Simulated Alternator Speed with 0.01 Sec. Dead Time 

Y. Simulated Engine Torque with 0.01 Sec. Dead Time 

Z. Simulated Alternator Speed - Loop Gain = 1.0 Per Unit 

AA. Simulated Alternator Speed - Loop Gain = 1.6 Per Unit 

BB. Simulated Alternator Speed - Loop Gain = 2.3 Per Unit 



ED .h ‘ - 



! 
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DISCUSSIQN OF RESULTS 

ALTERNATOR SPEED A comparison of strips A and B shows that there is 
excellent qualitative correlation between the simulation and the prototype 
performance. In the simulation the speed transient which occurs when the 
induction motor switches to speed five occurs approximately one second too 
late. This is due to an error in setting P90 which controls the switching 
time for motor speed five and in no way affects the validity of the model. 
The three criteria listed below give an indication of the quantitative cor- 
relation for the first speed 'transient. 

(1) Maximum deviation in speed — For the prototype the drop in fre- 
quency is eighteen chart divisions at a sensitivity of 20 mv/division. The 
frequency converter sensitivity is 2.5V/5Hz. The frequency deviation is 
0.72 Hz which corresponds to a 0.24 rev/sec drop in alternator speed. In 
the simulation the speed drop was 1.5 rad/sec or 0.239 rev/sec. 

(2) Ratio of speed drop to speed overshoot — Comparing the simulation 
and prototype gives the following approximate values : 

Simulation - 3.34 
Prototype - 3.60 

(3) Time to return to original speed--For both simulation and proto- 
type the speed first returned to synchronous speed in approximately 0.70 
seconds . 

The correlation noted above is admittedly a retrofit to the prototype 
data. However, it is clearly shown that the model used closely predicts 
the prototype performance. In the simulation the engine parameters set 
on the potentiometers were the values determined in the previous section. 
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The only settings which varied were the governor loop gain and the gain 
on the governor feedforward signal. These two settings are varied in the 
simulations to show their effects on engine performance . 

The simulated output of the feedforward network is shown on strip Q. 

As expected this signal shows the derivative operation in the FFWD net- 
work. Strips M through P show the effect of increasing the FFWD gain. 

On strip P the speed shows an initial increase when the load is applied. 
This type of response is not possible for the system under study so that 
an upper limit must be imposed on Q07 which corresponds to the FFWD gain. 
However, the trend in engine performance as Q07 is increased is evident. 

The effect of the FFWD on engine torque is shown in strips K and L. 

The effect of loop gain on performance is shown in strips Z, AA and 
BB. As gain is increased the speed deviations become smaller. In strips 
AA and BB it is evident that the FFWD gain is set too high since the speed 
shows a slight initial increase as load is applied. However, the influence 
of loop gain is still clear. 

In deriving the turbine model, the engine dead time was assumed to be 
negligible or at least less than 0.01 seconds. The influence of dead time 
on response is shown in strips U through Y. The dead time is simulated us- 
ing a first order lag approximation. In strip U the speed response for 0.1 
seconds dead time is shown. The response is noticeably worse than that 
shown on strip N which is the response for the same potentiometer settings 
for all components and with no dead time. The reason for the deterioration 
in response is evident from comparing strips V and W. The dead time causes 
the torque to lag the fuel actuator angle as shown. Strips X and Y corre- 
spond to a dead time of 0.01 second. Comparison of the speed response 



-78- 



shown on strips X and N shows that there is no noticeable difference in 
response. As a result, any dead time less than 0.01 second can be neglected. 

ALTERNATOR LOAD Strip C is the alternator KW load as measured by the 
governor load sensing section and as viewed through the 0.068 sec. LPF. 

Since the maximum alternator speed variation is on the order of one percent 
of base speed the load is roughly proportional to the alternator torque. 

The load signal was not calibrated to a known level of load so that this 
measurement may be used only for a qualitative comparison. Two ratios are 
used below to give an indication of the degree of correlation. The over- 
shoot ratio compares the maximum torque for each transient to the steady 
torque level at the corresponding induction motor speed. The steady state 
ratio compares the steady torque at each induction motor speed to the steady 
torque at the maximum motor speed. 





Table 


3. 








OVERSHOOT 


RATIO 


STEADY 


RATIO 


IM SPEED 


SIM 


PROTO 


SIM 


PROTO 


1 


2.80 


2.00 


.47 


.55 


2 


1.29 


1.24 


.65 


.70 


3 


1.16 


1.14 


.72 


.78 


4 • 


1.25 


1.19 


.84 


.87 


5 


1.61 


1.37 


1.00 


1.00 



The differences between simulated and observed values is due primarily to the 
difference between estimated and actual rotor resistance values. 

A rough check can be made on the actual magnitude of the simulated 
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alternator torque. At speed five the induction motor torque is approximate- 
ly 1400 lb-ft at 1150 RPM. Since the alternator speed is approximately 1200 
RPM and the power levels of the alternator and induction motor are approxi- 
mately equal, the alternator torque will be approximately 1400 lb-ft. The 
scaled value for the simulated alternator torque is 1350 lb-ft. 

INDUCTION MOTOR The correlation of the induction motor performance 
and the simulation is given in Table 4 below. The values for the motor 
torque and current were taken from the motor performance curve, Figure 11 . 
Also shown on this curve are the simulated values for torque and current. 

Table 4. 



Speed 

Position 


RPM 
Obs . 


RPM 

Sim. 


Torque 

(Performance 


1 RMS 
Curve ) 


Torque ■'■RMS 

Simulated 


i 


834 


840 


800 


- 


700 


200 


2 


966 


960 


1000 


310 


915 


250 


3 


10 32 


1015 


1100 


340 


1030 


275 


4 


1104 


1080 


1270 


380 


1190 


300 


5 


1180 


1150 


1450 


425 


1380 


350 



As with the engine simulation the results can vary widely depending on the 
values used for the motor and propeller parameters. The original estimates 
for the motor inertia and the rotor resistances were changed slightly in 
order to provide a better fit between the model and the observed performance. 
The differences between the original parameter estimates and the values used 
in the simulation are summarized in Table 5. 
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Table 5. 



Parameter 


Estimated 

Value 


Potentiometer 


Setting From 
Scaled Eqn. 


Setting Used 


^sm 


.71 ohms 


P10 


8460(10 3 ) 


8277(10 3 ) 


R a 


.0117 


P03 


0220 


0597 


R a(l) 


.240 


POO 


4450 


6000 


R a(2) 


.088 


P32 


1660 


2000 


R a(3) 


.072 


P65 


1360 


2000 


R a(4) 


.063 


P95 


1190 


2000 


Im + Ip 


21.6 


P35 


1475(10) 


2000(10) 




1.08 x 10~ 3 


P67 


3880 


3880 



All parameters agree with the estimated values except for the rotor 
resistances and the motor inertia. The difference in inertia is not large 
and is within the accuracy of the original estimate. The potentiometer 
settings corresponding to values of R a were considerably higher than the 
original estimates. These higher settings result from the fact that a 
single-phase model is being used to represent a three-phase machine. 

For the three phase delta connected rotor the resistance is given as 
0.0117 ohms. This resistance was probably measured across the slip rings 
and would correspond to the following diagram. 




As a result the single phase equivalent -resistance (Rp) is 1.5 R^.. This 
requires that the settings for the potentiometers representing the rotor 
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resistances in the table above be increased by a factor of 1.5. The values 
given in the manufacturer’s drawings for the switched resistors are assumed 
to be nominal values. For this simulation the same value (.70 ohms) was 
used for resistors two, three and four with resistor one assumed to be three 
times as large. These values are very close to the specified nominal values 
and give excellent qualitative correlation for the transient currents. 

With for the rotor given as 0.0117 ohms and corrected to a single 
phase equivalent the required setting on P03 is 0032 whereas a setting of 
0597 was used in the simulation. A value of 0032 was used but the current 
surge when switching to speed five was larger than observed on the prototype. 
As would be expected, the current transients at the other switching points 
were not noticeably affected by this change in the residual rotor resistance. 
The difference between the required rotor resistance potentiometer setting 
and the value actually used is due to the fact that at this point the divider 
A03 which is generating the motor back emf is poorly scaled for Analog compu- 
tation and accuracy is limited. 

When P65 is switched out of the motor resistance network the operation 
at A0 3 is as follows: 

Ra I 

v = 

N s - N 

At the switching point (N s - N) is approximately 0.12 per unit and approaches 
0.05 P.U. in the limit. With (I) properly scaled to avoid saturation at all 
other transients the magnitude of (I) at switching is approximately 0.3 P.U. 
With R a set at 0300 the component A03 is dividing a small signal by a small 
signal. As a result accuracy at this point is limited. A03 can be rescaled 
to give higher accuracy when N approaches N s but this will cause saturation 
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when the motor is first started. The divider must be constrained to 
operate at all times with the numerator smaller in magnitude than the 
denominator. 

Even with the reduced accuracy, when the motor speed approaches 
synchronous speed, there is still good correlation between the simulation 
and the prototype performance throughout the entire operating range. 

Strips E through H show the correlation of speed and current. The data 
tabulated above is plotted on the motor performance curve. Figure 11, and 
shows the quantitative correlation for torque and current. The simulated 
current shows the same trend as the performance curve data but is lower at 
all points. This is due in part to the error in measuring the magnitudes 
of the current from the strip recording. 
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7 . CONCLUSIONS AND RECOMMENDATIONS 

The correlation between the simulation and prototype for all observed 
variables proves the validity of the models used in the simulation. The 
natural tendency is to expect that highly complex models will be required 
for each component. However, reasonable accuracy in quantitative predic- 
tions and highly valuable qualitative information about system performance 
have been obtained using simple "off the shelf" models. The simulation 
can be helpful in design or as a training aid. 

As indicated earlier, the simulation gives the designer an indication 
of trends in system performance as the various parameters are adjusted. For 
a given motor the parameters X sm and R a are fixed. The switched resistors 
and switching times are the two most important parameters which can also be 
easily adjusted. The simulation shows how the performance can be varied 
over a wide range by adjusting these parameters. Although the designer has 
little control over the inertia of the motor and propeller, the inertia ef- 
fects are clearly shown. Finally, the effects of propeller pitch as char- 
acterized by the propeller torque coefficient can be investigated. 

The simulation provides an excellent study of the governing system 
and engine performance. The most important governor adjustments are con- 
trollable by potentiometer settings so that changes in the governor coef- 
ficients can be investigated. The effects on engine torque and speed 
response are clearly shown. In addition to providing useful design informa- 
tion, the simulation can demonstrate to operating personnel the effects on 
system performance due to changes in the governor settings. 

The transfer function approach to modeling the prime mover yielded a 
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simulation which can easily be adapted to any other type of prime mover by . 
substituting the proper transfer function (on a diode function generator 
if necessary) in the place of the gas turbine model. The model used in 
this simulation is readily adaptable for use with diesel engines with the 
most significant change being in the assumed engine dead time. 

The two most important dynamic effects which were not included in the 
simulation because of the limitation on computing hardware were the voltage 
regulation and fuel control valve effects. With additional hardward these 
effects can easily be included. However, as the results indicate, neglect 
of these effects did not seriously affect the accuracy of the simulation. 

Neglecting the stabilizer section of the governor control box and 
dynamics of the fuel control valve did not seriously affect the results 
for two reasons. The input disturbances to the gas turbine were below the 
frequency level where the stabilizer section influences performance. Al- 
though the load transient was close to a step change the second order filter 
on load inputs to the gas turbine limits change to approximately 12 Hz. Al- 
so, in scaling the simulation the developed engine torque was limited to 
3000 lb-ft. This provides a form of acceleration limiting. 

Neglect of the voltage regulation effects can influence the results 
in two ways. First, in the simulation the induction motor is supplied with 
a constant voltage whereas in the prototype the voltage varies as shown in 
Strip S. The error caused by this assumption is not serious. The second 
effect of neglecting the voltage regulation is on alternator torque. With 
field current assumed constant, the gyrator modulus on the alternator is 
constant so that the torque is directly 'proportional to the line current. 
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The field current does increase slightly during the transients so that the 
prototype alternator torque is higher than shown on strip J. However, the 
correlation between the simulation and prototype performance shows that 
considering field current constant does not cause a serious error. 

The use of a single-phase model for the three-phase machine intro- 
duces no serious errors. The inertia of the induction motor is large 
enough so that the double frequency torque pulsations in the motor are not 
noticed. At the alternator end, the second order filter provided adequate 
filtering of the pulsations and adequate frequency response in order to pass 
the transient loads. Even with unlimited Analog hardware it would not be 
necessary to use a three-phase model as the single-phase model is more than 
adequate . 



RECOMMENDATIONS 



The observed system performance indicates that the original system 
design is more than adequate in order to meet its performance specifica- 
tions. Therefore, use of this simulation in an attempt to improve the 
overall design is not practical. The only change to the original system 
which could be made at small cost is in the starting resistors. The speed 
range of the motor with the installed resistors is from approximately 65 
to 95 percent of synchronous speed. This range could be extended by chang- 
ing the resistors. The required resistors and dynamic performance can be 
determined by the simulation. However, it is not practical to lower the 
speed significantly because then the motor will be operating at a very low 



efficiency . 
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One area where additional work may be done is in utilizing the hybrid 
capabilities of the computer facility. Although the simulation has been 
done with only the Analog computer, the hybrid facility can be utilized to 
increase the effectiveness of the simulation. One very interesting study 
would be to perform an optimization study of the system. Even though 
prototype system performance is adequate, such a study would give a student 
the opportunity to learn the basics of hybrid computing and to study in the 
area of optimal control. 

The most important dynamic elements which were neglected were the fuel 
control valve and the voltage regulator. The governor control box stabilizer 
section was also neglected. Since the simulation required approximately 90 . 
percent of the Analog hardware, it would be difficult to include reliable 
models of these components. A more detailed study of the generator set 
alone can be performed by replacing the induction motor load by a step load 
change. It will then be possible to include the effects of the fuel control 
valve, voltage regulator and stabilizing section. 

The generator set frequency regulation is extremely important for 
proper operation of the shipboard electronic equipment. As a result, the 
speed deviations of the gas turbine must be kept to a minimum when the 
generator is supplying power to the electronic equipment. A simulation of 
the generator set in response to step load changes can be used to design 
the governing system and to determine the governor settings required for 
best performance. 

The most direct benefit of the simulation would be to use it to demon- 
strate to operating personnel how the system performance can be altered by 
adjusting the governor settings. The desire to "leave well enough alone" 
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together with a lack of complete understanding of how a system operates 
will usually prevent operators from adjusting the governor to get better 
performance. The simulation can be used to instruct operators in how the 
various governor settings effect performance. With a better understanding 
of how the system behaves, it will be easier to adjust the governor set- 
tings for near optimum performance. 
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